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Abstract 

The elastic and inelastic high-energy small-angle electron-positron scattering is con- 
sidered. All radiative corrections to the cross-section with the relative accuracy 5a jo = 
0.1% are explicitly taken into account. According to the generalized eikonal representa- 
tion for the elastic amplitude, in higher orders only diagrams with one exchanged photon 
may be considered. Single photon emission with radiative corrections and next-to- 
leading two-photon and pair production diagrams are evaluated, together with leading 
three-loop corrections. All contributions have been calculated analytically. We integrate 
the calculated distributions over typical for LEP 1 experiments intervals of angles and 
energies. To the leading approximation, the results are shown to be described in terms 
of kernels of electron structure functions. Some numerical results are presented. 
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1 Introduction 



An accurate verification of the Standard Model is one of the primary aims of LEP [|TJ. While 
electroweak radiative corrections to the s-channel annihilation process and to large-angle 
Bhabha scattering allow a direct extraction of the Standard Model parameters, small an- 
gle Bhabha cross-section affects, as an overall normalization condition, all observable cross- 
sections and represents an equally unavoidable condition toward a precise determination of 
the Standard Model parameters. The small-angle Bhabha scattering process is used to mea- 
sure the luminosity of electron-positron colliders. At LEP an experimental accuracy on the 
luminosity of 

| — | < 0.001 (1) 
(J 

has been reached 0. However, to obtain the total accuracy, a systematic theoretical error 
must also be added. This precision calls for an equally accurate theoretical expression for the 
Bhabha scattering cross-section in order to extract the Standard Model parameters from the 
observed distributions. An accurate determination of the small-angle Bhabha cross-section 
and of the luminosity directly affects the determination of absolute cross-sections such as, 
for example, the determination of the invisible width and of the number of massless neutrino 
species N u ||. 

In recent years a considerable attention has been devoted to the Bhabha process 0, [5|, |(| . 
The reached accuracy is however still inadequate 0. According to these evaluations the 
theoretical estimates are still incomplete; moreover, they are somewhat larger (~ a factor 
2) than the projected theoretical and experimental precision || and are comparable to the 
currently published experimental precision. 

The process that will be considered in this work is that of Bhabha scattering when electrons 
and positrons are emitted at small angles with respect to the initial electron and positron 
directions. We have examined the radiative processes inclusively accompanying the main 
e + e~— > e + e~ reaction at high energies, when both the scattered electron and positron are 
tagged within the counter aperture. 

We assume that the center-of-mass energies are within the range of the LEP collider 2e = 
y/s = 90 - 200 GeV and the scattering angles are within the range 6 ~ 10 - 150 mrad. We 
assume that the charged-particle detectors have the following polar angle cuts: 

B x <B-=p^q x =B<h, 9 2 < 9 + = P tf 2 < 6 4 , 0.01 < t < 0.1 rad , (2) 

where p 1; q x (p 2 , q 2 ) are the momenta of the initial and of the scattered electron (positron) 
in the center-of-mass frame. 

In this paper we present the results of our calculations of the electron-positron scattering 
cross-section with an accuracy of O(0.1%). The squared matrix elements of the various 
exclusive processes inclusively contributing to the e + e~ — > e + e~ reaction are integrated in 
order to define an experimentally measurable cross-section according to suitable restrictions 
on the angles and energies of the detected particles. The different contributions to the electron 
and positron distributions, needed for the required accuracy, are presented using analytical 
expressions. 

In order to define the angular range of interest and the implications on the required ac- 
curacy, let us first briefly discuss, in a general way, the angle-dependent corrections to the 



We consider e + e~ scattering at angles as defined in Eq. Within this region, if one 
expresses the cross-section by means of a series expansion in terms of angles, the main con- 
tribution to the cross-section da/d8 2 comes from the diagrams for the scattering amplitudes 
containing one exchanged photon in the t-channel. These diagrams, as well known, show a 
singularity of the type Q~ A for 9 — > 0, e.g. 

^ ~ 9^ . 

de 2 

Let us now estimate the correction of order 6 2 to this contribution. If 

^ ~ 9-\l+ Cl 9 2 ) , (3) 
then, after integration over 9 2 in the angular range of Eq. (Ef), we obtain: 



f)2 
max . 

J Z*9 2 ~ 9-1(1 + Cl 9ljn^). (4) 

g2 min 
min 

We see that, for 9 min = 50 mrad and # max = 150 mrad (we have taken the case where the 
9 2 corrections are maximal), the relative contribution of the 9 2 terms is about 5 x 10 -3 ci. 
Therefore, the terms of relative order 9 2 must be kept only in the Born cross-section where 
the coefficient c\ is not small. In higher orders of the perturbative expansion the coefficient 
Ci contains at least one factor a /it and therefore these terms can be safely omitted. This 
implies that, within our accuracy, only radiative corrections from the scattering-type diagrams 
contribute |7| . Furthermore only diagrams with one photon exchanged in the t-channel should 
be taken into account according to the generalized eikonal representation see Eq.(16) below. 

Having as a final goal for the experimental cross-section the relative accuracy of Eq. ([I]), 
and taking into account that the minimal value of the squared momentum transfer Q 2 = 
2e 2 (l — cos 9) in the region defined in Eq. @ is of the order of 1 GeV 2 , we may omit in the 
following also the terms appearing in the radiative corrections of the type m 2 /Q 2 , with m 
equal to the electron mass. 

The contents of this paper can be outlined as follows. In Section 2 we discuss the Born 
cross-section da B by taking the Z' -boson exchange into account and we compute the correc- 
tions to it due to the virtual and real soft-photon emission. We define also an experimentally 
measurable cross-section cr oxp with the experimental cuts on angles and energies taken into 
account and we discuss how to obtain it from the differential distributions. We present the 
results, as discussed above, in the form of an expansion in terms of the scattering angle 9. 
The ratio £ = cr cxp /cro is introduced by normalizing cr cxp with respect to the cross-section 
do = Ana 2 /e 2 9\. In Section 3, by using a simplified, but still suitable to obtain the required 
accuracy, version of the differential cross-section for the small-angle scattering, we discuss 
the contribution to cr cxp from the single bremsstrahlung process. The details of the Sudakov 



2 



technique which we use to calculate the hard-photon emission are given in Appendix A. In 
Section 4 we find all corrections of 0(a 2 ) to a cxp caused by virtual and real photons emission 
as well as pair production. In Section 5 we consider the virtual and soft-photon emission ac- 
companying the single photon bremsstrahlung process. The details of this derivation are given 
in Appendices B and C. In Section 6 we consider the double hard-photon emission process in 
both the same-side and opposite-side cases. Details are given in Appendix D. In Section 7 we 
consider the hard pair production process in both the collinear and semi-collinear kinematical 
region. The details of this calculation are given in Appendix F. In Appendix D the expres- 
sions for the leading logarithmic approximation in terms of structure functions factorization 
are given. The details of the cancellation of the A-dependence are presented in Appendix E, 
A is a small auxiliary parameter (A<1) which separate the processes of soft and hard real 
photon emission (see eq. (15)). In Section 8 the expressions to leading logarithmic 0(a 3 ) for 
the e + e~ and e + e~7 radiative processes are obtained. In Section 9, finally, estimates of the 
neglected terms together with numerical results are presented. 

A less detailed derivation of these results has been reported elsewhere ||. 



2 Born cross— section and one-loop virtual and soft cor- 
rections 



The Born cross-section for Bhabha scattering within the Standard Model is well known 
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Here 9 W is the Weinberg angle. In the small-angle limit (c = 1 — # 2 /2 + # 4 /24 + . . .), expanding 
formula (151) leads to 
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weak / 5 



(6) 
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where e = y/s/2 is the electron or positron initial energy and the weak correction term <5 wcak , 
connected with diagrams with Z°-boson exchange, is given by the expression: 

Q2 Q4 

<U k = igli - j(g 2 v + 9 2 a )Re X + —(gt + gt + Qg 2 v g 2 a )\x\ 2 ■ (7) 

One can see from Eq. (|7]) that the contribution c^ of the weak correction 5 wcak into the 
coefficient c\ introduced in Eq. (D) 



2 , (g 2 v +g 2 a ) M z a2 (g 4 v +gj + 6g 2 g 2 a ) M z 

S 2 g v H j -pr- + max ™ To" - 1 - W 

1 z 1 z 



"1 ~ ^at; ' A T-l ~ max or) p2 



The contribution connected with Z°-boson exchange diagrams does not exceed 0.3% for typical 
energy and angle ranges. Radiative corrections to Z° — 7 interference contributions were 
considered in detail in papers by W. Beenakker and B. Pietrzyk M. They are not small and 
should be taken into account in an analysis of experimental data. We will not touch the 
subject in this publication. 

In the pure QED case one-loop radiative corrections to Bhabha cross-section were calcu- 
lated a long time ago 0. Taking into account a contribution from soft-photon emission with 
energy less than a given finite threshold As, we have here for the cross-section daq ED , in the 
one-loop approximation: 

j (i) j B 

aa QED aa QED n . c . x n / n x 
dc = ^ (1 + Kirt + dsoft), (9) 

where dag ED is the Born cross-section in the pure QED case (it is equal to da B with 
9a = 9v = 0) and 



sin 2 (6>, 2) 

^virt "I" ^soft = 2 

7T 



2 n A 

2 fl-l n ZL + 2 ln(cot-)Vn-^+ / — ln(l-x) 
V m 2 2 / Ae J x 
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a 1 r 7T 2 
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H In , 

9 6 



23 11 , 4e 2 l a 1 r - 2 



7T (3 + C 2 ) 



^- (2c 4 - 3c 3 - 15c) 



Q 

+ 2 (2c 4 - 3c 3 + 9c 2 + 3c + 21) ln 2 (sin -) 

6 6 
- 4{c 4 + c 2 - 2c) In 2 cos - - 4 (c 3 + 4c 2 + 5c + 6) In 2 (tan -) 

2 8 9 

+ - (11c 3 + 33c 2 + 21c + 111) ln(sin -) + 2 (c 3 - 3c 2 + 7c - 5) ln(cos -) 

O Zi Zi 

+ 2 (c 3 + 3c 2 + 3c + 9) 5 t - 2 (c 3 + 3c) (1 - c) 5 S 

The value 5 t (5 a ) is defined by contributions to the photon vacuum polarization function U(t) 
(II(s)) as follows: 

. . a ( 1 5\ 1 ,q >0 

« = ; ^j'-i + i<;^ (10) 
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where 



In 



Q 2 



Q 2 = -t = 2e 2 (l-c), 



and we took into account the leading part of the two-loop contribution in the polarization 
operator. In the Standard Model, 5t contains contributions of muons, tau-leptons, VT-bosons 
and hadrons: 



6t = 6? + 8l + 6? + 6?, 5 s = 5 t (Q 



H 



(12) 



the first three contributions are theoretically calculable and can be given as: 
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- Vw {Vw - 4 ) + _ 
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The contribution of hadrons cannot be calculated theoretically; instead, it can be given as 
integration of the experimentally measurable cross-section: 



Q 2 



+?° e+e- 
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(14) 



For numerical calculations we will use for H(t) the results of Ref. [|Kj . 

In the small scattering angle limit we can present (P) in the following form: 



ao QED 

dc 
6 = 2* 
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This representation gives us a possibility to verify explicitly that the terms of relative order 
9 2 in the radiative corrections are small. Taking into account that the large contribution 
proportional to In A disappears when we add the cross-section for the hard emission, we 
can verify again that such terms can be neglected. Therefore we will omit in higher orders 
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the annihilation diagrams and multiple-photon exchange diagrams in the scattering channel. 
The second simplification is justified by the generalized eikonal representation for small-angle 
scattering amplitudes. In particular, for the case of elastic processes we have [11|: 



-1 i<p(t) 



7T S 



s > Q 2 > m 2 , (16) 



where A (s,t) is the Born amplitude, Fi(t) is the Dirac form factor and tp(t) = —a \n(Q 2 /X 2 ) 
is the Coulomb phase, A is the photon mass auxiliary parameter. The eikonal representation is 
violated at a three-loop level, but, fortunately, the corresponding contribution to the Bhabha 
cross-section is small enough (~ a 5 ) and can be neglected for our purposes. We may consider 
the eikonal representation as correct within the required accuracy^. 

Let us now introduce the dimensionless quantity £ = Q\ a cxp /(47ro; 2 ), with Q\ = e 2 9 2 , 
where <r cxp is the Bhabha-process cross-section integrated over the typical experimental energy 
and angular ranges^: 

O 2 r r r r ^^e+e-^e+iq^^e-iq^x^+x 

s = h e <^ - *•> e ? J A2 i> e > ^^dv ■ (17) 

where x± t 2, q% 2 are the energy fractions and the transverse components of the momenta of the 
electron and positron in the final state, sx c is the experimental cut-off on their invariant mass 
squared and the functions 0^ do take into account the angular cuts (0): 

e; = e(^3 - ^) - e 1 ), e c 2 = ®{e 4 - ^) e(^ - e 2 ). (is) 

X\E X\E x 2 e x 2 e 

In the case of a symmetrical angular acceptance (we restrict ourselves further to this case 
only) we have: 

02 = 0i, e 4 = e 3} p = ^>i. (19) 

We will present £ as the sum of various contributions: 

£ = £ + £ 7 + £ 27 + £ e+e ~ + £ 37 + S e+e " 7 (20) 
= E 00 (l + 5 + P + 5 27 + 5 e+e ~ + 5 37 + 5 e+e " 7 ), 
£ 00 = 1 - p~ 2 , 



4 Result obtained in paper [[L2|, we believe, is incorrect. It contradicts to the well established result of 
D. Yennie et al. about cancelation of infrared singularities. 

5 Really this quantity corresponds to some ideal detectors. It is intended for comparisons with the results 
of Monte Carlo event generators. 
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where S stands for a modified Born contribution, E 7 for a contribution of one-photon emission 
(real and virtual) and so on. The values of the 5 l as function of x c are given in Table 1 (see 
in Section 9). Being stimulated by the representation in Eq. (16), we shall slightly modify 
the perturbation theory, using the full propagator for the t-channel photon, which takes into 
account the growth of the electric charge at small distances. By integrating Eq. (j^j with this 
convention, we obtain: 

Zo = efj ^r(i - n(t))- 2 + s w + Eg, (21) 

el 

where T*w is the correction due to the weak interaction: 

el 

f d9 

T, w = Q x J -^5 wcak , (22) 

and the term E# comes from the expansion of the Born cross-section in powers of 6 2 , 

p 2 

= ^jf y (1 " n(-,Q?))- 2 + z9l^) . (23) 
i 

The remaining contributions to E in (^) are considered below. 



3 Single hard-photon emission 



In order to calculate the contribution to E due to the hard-photon emission we start from the 
corresponding differential cross-section written in terms of energy fractions x± }2 and transverse 
components q± 2 of the final particle momenta [ H |: 



where 
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(25) 
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and we use the full photon propagator for the t-channel photon. Performing a simple azimuthal 
angle integration of Eq. (24) we obtain for the hard-photon emission the contribution E^: 



l-A 2 

- / ±x-^F{x,D u D z -D 2 ,D,) 

TT J 1 — X 



(26) 



with 
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ai = (z 1 -z 2 ) +4z 2 a, a 2 = (zi - x z 2 ) + Ax z 2 a , a = — j(l — a; 



and the integration limits in (27) in the symmetrical case are: 

D x = x 2 , D 2 = 1, ^3 = x 2 p 2 , D 4 = p 2 . 
From Eqs. (^)-(|29|) we have that: 



(27) 



(2f 



(29) 



l-A 



2\\-2 



k(x, z) 



= - / dx- / — l-n {-zQl)) 

TT J 1 — X J Z 

x c 1 

x {[l + e(x 2 p 2 -z))(L-l) + k{x,z)}, 

- ' ' ~ [1 + Q(x 2 p 2 - z)] + U + 0(x 2 p 2 - L 2 + 9(z - x 2 p 2 )L 



(30) 



1 + x 2 

where L = ln(zQl/m 2 ) and 

x 2 (z — l)(p 2 — z 



Li 
L 3 



In 
In 



(x — z)(xp 2 — z) 
(z — x 2 ){xp 2 — z) 



111 



[z — x 2 )(x 2 p 2 — z) 



x 2 (x — z)(xp 2 — z) 



(31) 



(x — z)(x 2 p 2 — z) 



It is seen from Eq. ( |30"D that S H contains the auxiliary parameter A. This parameter disap- 
pears, as it should, in the sum S 7 = H H + TX +S , where Y, v+S is the contribution of virtual 
and soft real photons which can be obtained using Eq. (15): 



S7 = ^J^I d ^-^(-^Qi)r 2 {(L-i)p(x) 



(32) 



1 X c 



x [1 + e(x 2 P 2 - Z)\ + i±£ijfe(x, Z) - 5(1 - x) 

1 — X 
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where 

P(x) = (^) + = Inn {y^ 5(1 - x - A) + (2 + 21n A) 5(1 - x)} (33) 
is the non-singlet splitting kernel (see Appendix A for details). 

4 Radiative corrections to 0(a 2 ) 

A systematic treatment of all 0(a 2 ) contributions is absent up to now. This is mainly due 
to the extreme complexity of the analysis (more then 100 Feynman diagrams are to be taken 
into account considering elastic and inelastic processes). Nevertheless in the case of small 
scattering angles we may restrict ourselves by considering only diagrams of the scattering 
type. It is enough to make some rough estimates of other contributions. Contributions of pure 
annihilation-type diagrams, describing some 0(a 2 ) RC, have so-called double-logarithmical 



enhancement |TI| but, fortunately, it is proportional to the 4th power of the small scattering 
angle: 



{^)annih ~ (^lannih ~ 6\-f C\ (34) 

TT 

The contribution of interference of the scattering-type and the annihilation-type amplitudes 
can be estimated as 

(S^) mter/ ~ (£ e+e ") inter/ ~ # 2 (^) 2 ln 4 (^). (35) 

The absence of large lograrithms here has a similar origin as in (15). The uncertainties comes 
from the discussed contributions are considered in sect. 9. 

We consider first virtual two-loop corrections da vv to the elastic scattering cross-section. 
Using the representation (Jl6|) and the loop expansion for the Dirac form factor Fi 

F 1 = 1 + -F^ + (-) 2 if > (36) 

7T 7T 



one obtains 



^vv _ too&tfx _ n(t))- 2 [6(F 1 (1) ) 2 + 4if } ]. (37) 



dc dc TT 

The one-loop contribution to the form factor is well known: 



FP» = (L-1)1„A + 3 L _1 L2 _ 1+ 1 C2 . (3 , 
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The two-loop correction can be obtained from the results of Ref. 
form 

F (2) = tfrr + Ff e ~, 

where the contribution F± +e is related to the vacuum polarization by e + e~ pairs: 

19 
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The photon mass A entering Eqs. (|38|)-(Jl) is cancelled in the expression da^/dc for the 
sum of the virtual and soft-photon corrections of the second order d<jy V /dc (see Eq. fl3~7|)), 
dcr^g/dc and da^y/dc. 

The cross-section da S g/dc for the emission of two soft photons, each of energy smaller 
than Ae = eA, is (A < 1): 



da 
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7T 
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(42) 



and the virtual correction da sv /dc to the cross-section of the single soft-photon emission is: 



da 



fl = da A 2 (l-n(t))- 2 16if> 

7T 
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4 2 s 



(43) 



The contribution to X of this sum, except the part coming from Ff e connected with the 
vacuum polarization, contains no more than a second power of L. It has the following form: 



V77 
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It is convenient to separate the Rg+y m the following way: 



p77 
n S+V 

' s+v 



„77 



S+V + r S+V- T y + r 5+y 7 



rs+y 77 = L I 2 In A + 3 In A + 



9 



45 \ 



+ L — 4 In A — 7 In A + 3( 3 C2 
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(45) 
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The contribution to S coming from F^ +e contains an L 3 term, which is also cancelled 
when we take into account the soft pair production contribution 



71 



da (1 



2 r>e+e~ 
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-(L + 21nA) 2 
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(L + 21nA) 3 (46) 



27 



C 2 (L + 21nA)+0(l) 



Thus for the contribution of the virtual and soft e + e pairs to X we have 



^S+V 



pe+e 
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E? e e 



+ 4Ff 
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2D 4 
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In expressions (^5|) , A = 5e/e is the energy fraction carried by the soft pair, and it is 
assumed that 2m <C Se <C e. Here we have taken into account only e + e~ pair production. The 
order of magnitude of the radiative correction due to pair production is less than 0.5%. A rough 
estimate of the muon pair contribution gives less than 0.05% since ln(Q 2 /m 2 ) ~ 3 ln(Q 2 /m 2 ). 
Contributions of pion and tau lepton pairs give corrections that are still smaller. Therefore, 
within the 0.1% accuracy, we may omit any pair production contribution except the e + e~ 
one. 



5 Virtual and soft corrections to the hard-photon emis- 
sion 

By evaluating the corrections arising from the emission of virtual and real soft photons which 
accompany a single hard-photon we will consider two cases. The first case corresponds to the 
emission of the photons by the same fermion. The second one occurs when the hard-photon 
is emitted by another fermion: 



da 



da H(s+v) + daH{s+v) + da H +v) + da £ 



H(S+V) 

In the case when both fermions emit, one finds that: 



(S+V) 
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2S 



H 
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TT 



(L-l)lnA+-L- 1 



(49) 
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where T, H is given in Eq. (|30"1). A more complex expression arises when the radiative corrections 
are applied to the same fermion line. Here the cross-section may be expressed in terms of the 
Compton tensor with an off-shell photon |L7|, which describes the process 



7 *(g) + e-(pi) -> e-(gi) + 7 (A;) + ( 7<oft ). 
In the limit of small-angle photon emission we have: 
da H(s+v) _ a 4 dxd 2 qid 2 qi 

T 



■[(B 11 (s 1 ,t 1 )+x 2 B 11 (t 1 ,s 1 ))h + T} : 



T n {s\M) +x 2 Tn(ti,si) +x(T l2 (s 1 ,t 1 ) + T 12 (t 1 ,si)), 



(50) 



(51) 



h 



2 L-ln 



-Ml 



- 1 )(21nA-lnx) + 3L-ln 2 x 
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where A = (Ae/e) <C 1, Ae: is the maximal energy of the soft photon, escaping the detectors, 
Bn(si,ti) = (— 4(<7^) 2 )/(siti) — 8m 2 /si is the Born Compton tensor component, and the 
invariants are: Si = 2qik, ti = —2pik, Ui = (pi — qi) 2 , Si + t\ + U\ — q 2 . 
The final result (see Appendix C for details) has the form: 



^{S+V) 



■<H{S+V) 



I 71 



p2 dz 1 r A dx(i + x 2 ) r f 1 . , 3 

' — L \ [ 2 In A - lnx + - 

1 — x 



(52) 



x [(L — 1)(1 + 0) + k(x, z)] + - In 2 x + (1 + 0)[— 2 + lnx — 2 In A] 



-L In x + 2 In A In x — In x ln(l — x) 



U 2 {x) 



In 2 x — Li 2 (l — x) 
-/find-.), 



x(l — x) + 4xlnx 
2(1 + x 2 ) 



(1-x) 2 
2(1 + x 2 ) 



where /c(x, z) is given in Eq. (|30|) and 6 = 6(x 2 p 2 — 



6 Double hard-photon bremsstrahlung 

We now consider the contribution given by the process of emission of two hard photons. We 
will distinguish two ) the double simultaneous bremsstrahlung in opposite directions 

along electron and positron momenta, and b) the double bremsstrahlung in the same direction 
along electron or positron momentum. The differential cross-section in the first case can be 
obtained by using the factorization property of cross-sections within the impact parameter 
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representation [18|. It takes the following form 



(see Appendix A): 



dxid qidx 2 d q 2 



7T" 



7T 



(1 _ n(-(k^Y))-'R(x i; qf, k x )R{x 2 - q$, -k x ), (53) 



where R(x;q ,k ) is given by Eq. fl25|) . The calculation of the corresponding contribution 
to £ is analogous to the case of the single hard-photon emission and the result has the 
form: 



iff 

J H 



oo 

1 ,a, 2 f dz 
~ 1 -rl 



4 7T' 







2" 



n(-ztf)) 



1-A 1-A 

2u_J ' dxi J dx 2 

Xc/xi 



l + x\ l + x\ 



Xi 1 - x 2 



^(x 1 ,z)^(x 2 ,z), (54) 



where (see Eq. (|31|) ): 

= (I-l)[9(z-l)9(/-2) + e(z-x 2 )e(pV-z)] 



(55) 



+ L 3 [-e(x -z) + e( 2 -iV)] + (L 2 + 

+ (L 1 + tL^-)e(2-l)0(p 2 - :) 



X) 



1 + x 2 



e(z-x 2 )Q(x 2 p 2 -z) 



1 + x 

+ (0(1- z)-e(z- p 2 ))ln 



(2; — x)(p 2 — z) 



(xp 2 — z) (z — 1) 



Let us now turn to the double hard-photon emission in the same direction and the hard 
e + e~ pair production. Here we use the method developed by one of us fl9|, |2(|. We will 
distinguish the collinear and semi-collinear kinematics of final particles. In the first case all 
produced particles move in the cones within the polar angles 61 < 9 <^ 1 centered along the 
charged-particle momenta (final or initial). In the semi-collinear region only one produced 
particle moves inside those cones, while the other moves outside them. For the totally inclusive 
cross-section, such a distinction no longer has physical meaning and the dependence on the 
auxiliary parameter 8q disappears. We underline that in this way all double and single- 
logarithmical contributions may be extracted rigorously. The contribution of the region when 
both the photons move outside the small cones does not contain any large logarithm L. The 
systematic omission of those contributions in the double bremsstrahlung and pair production 
processes is the source of uncertainties of order (a/n) 2 < 0.6 • 1CT 5 . 

The contribution of both collinear and semi-collinear regions (we consider for definiteness 
the emission of both hard photons along the electron, since the contribution of the emission 
along the positron is the same) has the form (see Appendix B for details): 



J HH — TV - 



<y ]2 

4 l 7T 



Z 2 



u(-zQl)y 



(56) 
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tHH 



1-2A l-is-A 

dx J dxi 

A 

„2 "2 



I HH L 



X\{1 — X — X\)(l — Xi)' 



A 0(xV - z) + B + C 0((1 - xi)V - 2 



where 
A = 

B = 



7/? (§ + In 



(p 2 x 2 — z)' 



x 2 (p 2 x(l — Xi) — z)' 



+ (x 2 + (l-xi) 4 )ln 



(1 — Xi) 2 (l — X — Xi) 



XX i 



x 2 - l)(p 2 - z)(z - x 2 )(> - (1 - xi) 2 ) 2 (p 2 x(l - Xi) - zf 



(p 2 x 2 — z)(z — (1 — xi)) 2 (p 2 (l — xi) 2 — z) 2 (z — x(l — Xi)) 2 



c 



where 



(x 2 + (l-x 1 ) 4 )ln ( * Xl)2x \ +S B , 

X{1 — X — Xi) 

x(p 2 (l — x x ) 2 — z) 2 



' ' ' 1 ' 2 |U (1 - xx) 2 (p 2 x(l - Xl ) - z)(p 2 (l - Xl ) - z) 
2(l-xi) / 3-2x(l -xi)7. 



7 
7a 



l + (l-xx) 2 , /3 = x 2 + (l-X!) 2 , 

xxi(l — x — Xi) — x 2 (l — x — Xi) 2 — 2(1 — Xi)P, 

xxi(l — x — Xi) — x 2 (l — x — Xi) 2 — 2x(l — Xi)7. 



One may see that the combinations 

r77 + yH(S+V) + yHH 



r l i yH , yA'+V , yli 



5+y , v# 



with r 77 and r 7 normalized (see Eqs. (42,43)) to 



+ 7A, 



(57) 



nil 



-? 'f%{l~Tl{-zQ\))- 2 r 

TV J Z 



and 



-) 2 /^(i-n(-zQ?))- a r2 +V7) 

71 J Z 



respectively, do not depend on A for A — > (see Appendix E). 

The total expression S 27 , which describes the contribution to ( f20l) from the two-photon 
(real and virtual) emission processes is determined by expressions (43), (47), (49), (51) , (53) 
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and (55). Furthermore it does not depend on the auxiliary parameter A and has the form: 

£ 2 ^ = Y%L V + 2S^ +S ) + 2Sf +y + Eg + 2E HH (5* 



e 2 e\ 



£77 + S 7 + (_)^(077 + 07) ; 



C = In 



77T 



The leading contributions £ 77 , S 7 have the following forms (see Appendix D): 



■'77 



2 71' 



J ^L\l - U(-Qlz))- 2 J dx { V 2 )(x) [ 6(xV - z) + 1] 

1 x c 

+ /jP(t)P(f)e(*V-^)}, 



(59) 



p(2), 



/ ^p( t ) p/£) = i im 

a; 

1 +X 2 



;r 



21nA + - 

2 



4C 2 



5(1 -x) 



(60) 



V7 



iQf J^L 2 (1-U(-Qlz))- 2 Jd Xl J dx 2 P( Xl )P(x 2 ) 



1 — x 

oo 



2 ln(l — x) — lna; H — I H — (1 + x) hxx — 1 



6(1 - x - A) j, 



dz 



(61) 



Xo/a;i 



x [0(^ - l)0(p 2 - z) + 9(z - x 2 )0(x 2 p 2 - z)] 
x [6(z - l)6(p 2 - z) + 9(z - a$6(a^p 2 - z)]. 

We see that the leading contributions to S 27 may be expressed in terms of kernels for the 
evolution equation for structure functions. 

The functions 77 and 7 in expression Eq. ( |5Tf ) collect the next-to-leading contributions 
which cannot be obtained by the structure functions method | 2l| . They have a form that 
can be obtained by comparing the results in the leading logarithmic approximation with the 
logarithmic ones given above. 



7 Pair production 

Pair production process in high-energy e + e~ collisions was considered about 60 years ago (see 



14| and references therein). In particular it was found that the total cross-section contains 
cubic terms in large logarithm L. These terms come from the kinematics when the scattered 
electron and positron move in narrow (with opening angles ~ m/e) cones and the created 
pair have the invariant mass of the order of m and moves preferably along either the electron 
beam direction or the positron one. According to the conditions of the LEP detectors, such a 
kinematics can be excluded. In the relevant kinematical region a parton-like description could 
be used giving L 2 and L-enhanced terms. 
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We accept the LEP 1 conventions whereby an event of the Bhabha process is defined as 
one in which the angles of the simultaneously registered particles hitting opposite detectors 
(see Eq. (||)). 

The method, developed by one of us (N.P.M.) JH5], pOfl , °f calculating the real hard pair 
production cross-section within logarithmic accuracy (see the discussion in sect. 6) consists 
in separating the contributions of the collinear and semi-collinear kinematical regions. In the 
first one (CK) we suggest that both electron and positron from the created pair go in the 
narrow cone around the direction of one charged particle [the projectile (scattered) electron 
p 1 (q-j) or the projectile (scattered) positron p 2 (q 2 )]'- 

P3>-~P3>i~P^Pi<0o<.l, 4/™ » 1, Pi =Pi, P 2 > 9i> <Z 2 (62) 

The contribution of the CK contains terms of order (aL/n) 2 , (a/7r) 2 Lln(9 /9) and (a/7r) 2 L, 
where 9 = pZq\ is the scattering angle. In the semi-collinear region only one of conditions 
(|6"2"D on the angles is fulfilled: 



p + P- < O , p±pi > 9 ; or pjp { < Q , p + p { > 9 ; (63) 
or p~ZPi > O , p^p i < 9 . 



The contribution of the SCK contains terms of the form: 

The auxiliary parameter 9$ drops out in the total sum of the CK and SCK contributions. 

All possible mechanisms for pair creation (singlet and non-singlet) and the identity of the 
particles in the final state are taken into account ||24|| . In the case of small-angle Bhabha 
scattering only a part of the total 36 tree-type Feynman diagrams are relevant, i.e. the 
scattering diagrams^. 

The sum of the contributions due to virtual pair emission (due to the vacuum polarization 
insertions in the virtual photon Green's function) and of those due to the real soft pair emission 
does not contain cubic (~ L 3 ) terms but depends on the auxiliary parameter A = 5e/e 
(m e C fc < e, where Se is the sum of the energies of the soft pair components). The A- 
dependence disappears in the total sum after the contributions due to real hard pair production 
are added. Before summing one has to integrate the hard pair contributions over the energy 
fractions of the pair components, as well as over those of the scattered electron and positron: 

A = < X\ + 22, x c < x = 1 — X\ — %2 < 1 — A, (65) 

X\ — , X2 — , X — — , 

see 



6 We have verified that the interference between the amplitudes describing the production of pairs moving 
in the electron direction and the positron one cancels. This is known as up-down (interference) cancellation 
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where e± are the energies of the positron and electron from the created pair. We consider for 
defmiteness the case when the created hard pair moves close to the direction of the initial (or 
scattered) electron. 

Consider first the collinear kinematics. There are four different CK regions, when the 
created pair goes in the direction of the incident (scattered) electron or positron. We will 
consider only two of them, corresponding to the initial and the final electron directions. For 
the case of pair emission parallel to the initial electron, it is useful to decompose the particle 
momenta into longitudinal and transverse components: 

p + = X\V\ + pt, P- = X2P1 +P-, qi = xpi + qi, (66) 
x = 1 - x 1 - x 2 , q 2 ~ P2, pi+p A l + qi = 0, 

where pf are the two-dimensional momenta of the final particles, which are transverse with 
respect to the initial electron beam direction. It is convenient to introduce dimensionless 
quantities for the relevant kinematical invariants: 

2 f n \ 2 

< Zi < — > 1, (67) 




A = = (xix 2 ) [(l-x) + x 1 x 2 (z 1 + z 2 - 2y/ziz 2 cos0)] 



in 



Ai = PlP — x 2 1 [1 + x\ + x\z 2 \ , A 2 = ^ PlP+ = x 1 1 [1 + x\ + x\zi\ 
m 1 m 2 



C 



{ J^^ = 2-A 1 , D= { J^^-l = A-A 1 -A 2 , 

where <fi is the azimuthal angle between the (P1P+) and (pxPZ.) planes. 

Keeping only the terms from the sum over spin states of the square of the absolute value of 
the matrix element, which give non-zero contributions to the cross-section in the limit #0 — > 0, 
we find that only 8 from the total 36 Feynman diagrams are essential p4 |. 

The result has the factorized form: 

El^l 2 „ =El^oi 2 2Wij, (68) 
p+,p- pi m 4 

spins spins 

where one of the multipliers corresponds to the matrix element in the Born approximation 
(without pair production): 

E IM0I 2 = 2W f g4 + f, 2 + " 4 ) , (69) 

spins \ / 

s = 2pip 2 , t = —Q 2 x, u — —s — t, 

and the quantity /, which stands for the collinear factor, coincides with the expression for I a 
obtained in [20]. We write it here in terms of our kinematical variables: 

, = (i-^( ^-^ + ^ y + (i-,)-'( c "-;),- c " )' (to) 
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+ 



+ 



+ 



2(i- X2 ) 2 -(i- x y 

■IrAD l-x 

i r(i - x 2 ) 2 - 2(1 - x) 2 



I XlX ~ X2 , Q I \ 

+ — V 3(x 2 - x) 



2xCD 



1 — x 2 



+ 



x 2 (x 2 + xl) 



2x(l - 
AC 



3x 

x)AC + D 2 + AD 2 



1 — x 2 

X — %-\ Xo 

-; 4 

1 — x 



+ 



AA 



1 



CD 2 
Xi - x)(l + x 2 



3(x 2 - 
2T1 



x(l - x 2 ) 



We rewrite the phase volume of the final particles as 

d 3 gid 3 g 2 



- 2- 



xA 2 D 
1 — x 



X 



dr 



X 



(27r)62g?2g 2 
m 4 2~ 8 



7T 



o(27r)V 4) (pix+p 2 -gi-g 2 ) 
d0 



Xix 2 dxidx 2 dzidz 2 — . 

27T 



(71) 



Using the table of integrals given in Appendix F we further integrate over the variables of the 
created pair. Following a similar procedure in the case when the pair moves in the direction 
of the scattered electron, integrating the resulting sum over the energy fractions of the pair 
components, and finally adding the contribution of the two remaining CK regions (when the 
pair goes in the positron directions), we obtain 7 : 



d. coU = ^/^L{^)(L + 21n^)(l + e) 
+ AR (x)\nx + 2Gf(x)+2f 1 (x)\, A = -J 

EE &(x 2 p 2 - Z) ~- 



J 1, x 2 p 2 > z, 
\ 0, x 2 p 2 < z, 

Ro(x) = ^I^ + ii— ^>(4 + 7x + 4x 2 ) + 2(l + x)lnx, 

107 136 2 2 4 20 2 r A 2 r 

tlx) = 1 x X ; H — — Ax — bx + 1 

JK ' 9 9 3 3x 9(1 -x) 3 l 

4 i i / \ 1 r *? 13 2 t O 

+ — lnfl — x) H — 8x + 5x — 7 lux In x 

x(l-x) J v ; 3 L l-x J l-x 

2(3x 2 - 1) 

+ 4(1 + x) In a; In (1 - x) -Li 2 (l-x), 

1 — x 

r / n „ f/ l, 116 127 4 2 2 20 2 r 2 

/i x = -xRe f (-) = -— + —x + -x 2 + — -— - + --4x 2 

x 9 9 3 3x 9 1 — x 3 



(72) 



7 Some misprints, which occur in the expressions for f(x) and f\(x) in pfl, p4j, are corrected here. 
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4 I 5 

- 5x + 1 + r] ln(l -x) + -\8x 2 - lOx — 10 H ] lnx 

x(l — x) 3 1 — a; 

2( x 2 - 3) 

- (1 + x)ln 2 x + 4(l + x)lnx\n(l - x) -Li 2 (l -x), 

1 — x 

Ql = £#min, L = 111 . 

Consider now semi-collinear kinematical regions. We will restrict ourselves again to the 
case in which the created pair goes close to the electron momentum (initial or final). A 
similar treatment applies in the CM system in the case in which the pair follows the positron 
momentum. There are three different semi-collinear regions, which contribute to the cross- 
section within the required accuracy. The first region includes the events for which the created 
pair has very small invariant mass: 

4m 2 < (p++p~) 2 < \q 2 \, 

and the pair escapes the narrow cones (defined by 6q) i n both the incident and scattered 
electron momentum directions. We will refer to this SCK region as p + \\ p . The reason 
is the smallness (in comparison with s) of the square of the four-momentum of the virtual 
photon converting to the pair . 

The second SCK region includes the events for which the invariant mass of the created 
positron and the scattered electron is small, 4m 2 <C (p + + qi) 2 <C \q 2 \, with the restriction 
that the positron should escape the narrow cone in the initial electron momentum direction. 
We refer to it as p + || q 1 

The third SCK region includes the events in which the created electron goes inside the 
narrow cone in the initial electron momentum direction, but the created positron does not. 
We refer to it as p_ || p 1 p4 . 



The differential cross-section takes the following form: 

a 4 \M\ 2 dx x dx 2 dx 2 ± 2 x 2 ± 2 ± 
da = — — —— dp+dp _ dqidqi 6(l- Xl -x 2 -x) (73) 



x 



where X\ (x 2 ), x and p^_ (pi), q\ are the energy fractions and the perpendicular momenta of 
the created positron (electron) and the scattered electron (positron) respectively; s = (P1+P2) 2 
and q 2 = —Q 2 = (p 2 — g 2 ) 2 = —e 2 9 2 are the center-of-mass energy squared and the momentum 
transferred squared; the matrix element squared \M\ 2 takes different forms according to the 
different SCK regions. 

For the differential cross-section in the p + || p_ region we have (see, for details, p3|): 

j a \ A A d^) 2 d(^) 2 

dap up = —L dx dx 2 74 
+ TT (qf) 2 (qf + qf) 2 



d(f) 

x 



2vr (qf + xq$ 



2 , 1 x 2 AXX1X2 

Xl) +{1 - X2) -Jl^x? 
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where <fi is the angle between the two-dimensional vectors and q 2 , qj; 2 are the transverse 
momentum components of the final electrons, Xi j2 are their energy fractions {x — 1 — X\ — x 2 ). 
At this stage it is necessary to use the restrictions on the two-dimensional momenta q\ and q 2 . 
They appear when the contribution of the CK region (which here represents the narrow cones 
with opening angle 8q in the momentum directions of both incident and scattered electrons) 
is excluded: 



Pa 



£+ 



>9 , 



P± 

£+ 



Si. 

£2 



(75) 



where e + and e 2 are the energies of the created positron and the scattered electron respectively. 
In order to exclude from the above equation we use the conservation of the perpendicular 
momentum, in this case: 



1 1 J- 1 

qi + Q2 + pi 

X\ 

In the semi-collinear region p + || q x we obtain: 

d(q^) 2 d(^) 2 



0. 



do 



PJQ, 



X 



— L dx dx2 

7T 

d(j) 



(qf 



(76) 



x 



2n (qi + xq^f (1 - x 2 f 



x 



Xij 



AxxiX 2 



(1-^2) 



with the restrictions 



£2 



0. 



x 



(77) 



Finally for the p_ || p x semi-collinear region we get: 

d(^) 2 d(^) 2 



do 



P Pi 



x 



— L dx dx 2 

7T 

d(j) 



IqfV (qty 

[1 - x) 2 + (1 — x x ) 



(7E 



2vr (qi + q^f 



AxxiX 2 



(1-^2) 



The restriction due to the exclusion of the collinear region when the created pair moves 
inside a narrow cone in the direction of the initial electron has the form 



\P+ 

£1 



>9o, 



pi + qf + q2= 0. 



(79) 



In order to obtain the finite expression for the cross-section we have to add dop + \\p_ 



d<Tp + ||g i +dop np i to the contribution of the collinear kinematics region ( [72] ) and those due to 
the production of virtual and soft pairs. Taking into account the leading and next-to-leading 
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terms we can write the full hard pair contribution including also the pair emission along the 
positron direction, after the integration over x 2 as 



C 2 {x) 



R{x) 



2^f/^ / dxi [ L 2 (l + &)R(x) + C[QF 1 (x) + F 2 (x)] 



Fi(x) = d{x) + Ci(x), F 2 (x) = d(x) + C 2 (x) 
d(x) 

Ci(ar) 



1^(1^-)-?)= 



113 142 2 2 4 4 . , , 
— + —x - -x - — - - (1 + x) ln(l - x) 



+ 



2 1 + x 2 f, (x 2 p 2 - 



3 1 - x 
7 



In 



(xp 2 — z) 2 



- 3Li 2 (l - x) 



1 — x 



lnx + 



2(5x 2 - 6) . 2 

' ln 2 x + (3(x) In 



1 — x 



f (8x 2 + 3x - 9 
(x 2 p 2 — z) 2 

? ' 



X 



122 133 4 2 2 4. n . , ^ . 
— + —x + -x + — - -(1 + x) ln(l - x) 



+ 



2 1 + x 2 r 

3 1 -x 



In 



(^-x 2 )(p 2 -^-l) 



(x 2 p 2 — z) (z — x) 2 

1 13 S 

+ -(-8x 2 - 32x - 20 



3Li 2 (l-x) 



H — ) In x + 3(1 + x) In 2 x 

1 — x x 



+ /3(x) In 



(z - x 2 )(p 2 - z)(z - 1) 



x 2 p 2 — 2 



= 2R(x) - 



2 1 + x 2 

3 1 -x 



\ \ + X + V-^( 4 + 7x + 4x 2 ) + (1 + x) lnx. 
3 1 — x ox 



*0) 



(81) 



Eq. ( p0|) describes the small-angle high-energy cross-section for the pair production pro- 
cess, provided that the created hard pair can move along both electron and positron beam 
directions. 

The contribution to the cross-section of the small-angle Bhabha scattering connected with 
the real soft (with energy lower than Ae) and virtual pair production can be defined ]23| by 
the formula: 



4a 4 

20 
"9 



dz 



L z [ =-lnA + I ) £ 



17 4 ,2a 

h - In A 

6 3 



(82) 



MA-- C2 )}. 

Using Eqs. (|80[) and (jS2"D it is easy to verify that the auxiliary parameter A is cancelled in the 
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sum cr pair = cr 



hard "T~ ^aoft+virt 



We can, therefore, write the total contribution cr pair as 



2a 4 r dz ( 0/ 4, , , 2 r dx -. 
^QlJ A L{1+ 3 Hl -^-3ll-^ e) + C 



nQ 2 
s 40 
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y 



-3C2 - — ln(l - O + - ln 2 (l - x c ) + J —0 • (— - - ln(l - x)) 

1 

+ J dx[L 2 (l + Q)R(x) + £(9(7x0) + C2O))]}, 5(a?) = R(x) - — - 



x 



(83) 



= 1-0. 

The right-hand side of Eq. (|83| ) gives the contribution to the small-angle Bhabha scattering 
cross-section for pair production. It is finite and can be used for numerical estimations. The 
leading term can be described by the electron structure function D e e [x) p2jj . 



8 Terms of 0{aCf 



In order to evaluate the leading logarithmic contribution represented by terms of the type 
(a£) 3 , we use the iteration up to /3 3 of the master equation J2j] obtained in Ref. [p2] . To 
simplify the analytical expressions we adopt here a realistic assumption about the smallness 
of the threshold for the detection of the hard subprocess energy and neglect terms of the order 
of: 



x n c (-Cf < 3 - 10" 

7T 



n = 1,2,3. 



(84) 



We may, therefore, limit ourselves to consider the emission by the initial electron and positron. 
Three photons (virtual and real) contribution to E have the form: 



1 P \] 1 I 

= i^ c ) 3 J^l dx iJ d ^ e ( x ^- x c) qS(1-x 2 ) P (3) 0n) (85) 



:i + o(xD), 



where P(x) and P^ 2 \x) are given by Eqs. (^) and (|59|) correspondingly: 



2 9 



0!0 2 = 0^-^) 0(p 2 ^-,), 

p( 3 )( x ) = 5(1 - x) A t + 0(1 - x - A) e t , 
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48 



48 



IC3-i C2 (lnA + ?) + i(lnA + f 3 



(86) 



1 l+i 2 

2 1 - x 
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1 



Co H — lnfl — x) Inx H — ln 2 (l — x) 

32 2 s 4 v ; 8 2 v ; 



+ — ln 2 x In x lnfl — x) 

12 2 K ' 



1 1 
+ -(1 + x) lnxln(l — x) — — (1 — x) ln(l — x) 



+ — (5 - 3x) Inx - — (1 - x) - ^(1 + x) In 2 x + ^(1 + x)Li 2 (l - x)|. 

The contribution to S of the process of pair production accompanied by photon emission 
when both, pair and photons, may be real and virtual has the form (with respect to paper by 
M. Skrzypek |22| we include also the non-singlet mechanism of pair production): 



1 ,a 



£) 3 dz z 2 dxi / dx 2 0(xix 2 — x c ) 



x{-[^(xO - -i? s (x!)] 5(1 - x 2 )9(x 2 p 2 - z) + ^ ^(^2)^1) OxG,}, 



where 



i?(x) = R s {x) + -P(x), i? s (x) 



1 — x 
3x 



(4 + 7x + 4x 2 ) + 2(1 + x) In x, (87) 



R p {x) = R s {x){- + 21n(l -x)) + (1 + x)(-ln 2 x + 4Li 2 (l - x) 



+ ~(-9 - 3x + 8x 2 ) Inx + -(— - 
3 3 x 



3x + 3x 2 ) + -P (2) (x) 
3 



The total expression for £ in Eq. (20) is the sum of the contributions in Eqs. (21), (32), (56), 
(60), (66) and (68). The quantity £ depends on the parameters x c ,p and Q\. 

9 Estimates of neglected terms and numerical results 

The uncertainty of our calculations is defined by neglected terms. Let us list them, 
a) Terms of the first order RC coming from annihilation-type diagrams (15): 



H*/£a.<o.io.i«- 



b) Similar terms in the second order do not exceed (see sect. 4) 



7T' 



9l J 6* 



39) 



7T 



2 ^ 4 -# 4 )£ 4 <0.5-10 



-5 
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c) We neglect terms which violate the eikonal approximation: 



a Q 2 K . . 

- — < 0.3 • 10~ 6 . 90 

71 S 

d) We omit term of the second order which are not enhanced by large logarithms: 

.a. 



^-0.5-10" 5 . (91) 
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e) Creation of heavy pairs ([ifi, tt, titt, . . .) gives in sum at least one order of magnitude 
smaller than the corresponding contribution due to light particle production pSfl : 



+ + S rT < 0.1 S e e ~ < 0.5 • 10~ 4 . (92) 

f) Higher-order corrections, including soft and collinear multi-photon contributions, can 
be neglected since they only give contributions of the type (aL/ir) 4 < 0.2 ■ 10~ 5 or less. 

g) The terms in the third order associated with the emission off the final particles^: 

x c {—) 3 < 0.3 ■ 10" 4 (for x c = 0.5). (93) 

71 

Regarding all the uncertainties a)-g) and (82) as independent ones we conclude the total 
theoretical uncertainty of our results to be ±0.006%. " 

Let us define S° to be equal to £o|n=o ( see Eq. (21)), which corresponds to the Born 
cross-section obtained by switching off the vacuum polarization contribution IT(t). For the 
experimentally observable cross-section we obtain: 

o = -^S° (1 + 5 Q + cF + 6 2 ~< + 5 e+e ~ + 5^ + 5 e+e ^), (94) 



Ql 



where 



and 



S° = E |n=o = l-p" 2 + SH/ + S,|n=o (95) 



y yO y-y y2-y 

Z^ Zjq Zj 



The numerical results are presented below in Table 1. 

Table 1: The values of 5 i in per cent for ^ = 91.161 GeV, 6 1 = 1.61°, 9 2 = 2.8°, 
sin 2 6 W = 0.2283, T z = 2.4857 GeV. 



^Usually, in a calorimetric experimental set-up such terms do not contribute. 
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x c 


So 




^leading 


r^ 7 

^nonlcading 


5 e+e ~ 




5 3 7 




0.1 


4.120 


-8.918 


0.657 


0.162 


-0.016 


-0.017 


-0.019 


-4.031±0.006 


0.2 


4.120 


-9.226 


0.636 


0.156 


-0.027 


-0.011 


-0.016 


-4.368±0.006 


0.3 


4.120 


-9.626 


0.615 


0.148 


-0.033 


-0.008 


-0.013 


-4.797±0.006 


0.4 


4.120 


-10.147 


0.586 


0.139 


-0.039 


-0.005 


-0.010 


-5.356±0.006 


0.5 


4.120 


-10.850 


0.539 


0.129 


-0.044 


-0.003 


-0.006 


-6.115±0.006 


0.6 


4.120 


-11.866 


0.437 


0.132 


-0.049 


-0.002 


-0.001 


-7.229±0.006 


0.7 


4.120 


-13.770 


0.379 


0.130 


-0.057 


-0.001 


0.005 


-9.194±0.006 


0.8 


4.120 


-17.423 


0.608 


0.089 


-0.069 


0.001 


0.013 


-12.661±0.006 


0.9 


4.120 


-25.269 


1.952 


-0.085 


-0.085 


0.005 


0.017 


-19.379±0.006 



Each of these contributions to cr has a sign that can change because of the interplay 
between real and virtual corrections. The cross-section corresponding to the Born diagrams 
for producing a real particle is always positive, whereas the sign of the radiative corrections 
depends on the order of perturbation theory. For the virtual corrections at odd orders it is 
negative, and at even orders it is positive. When the aperture of the counters is small the 
compensation between real and virtual corrections is not complete. In the limiting case of 
small aperture (p — > 1, x c — > 1) the virtual contributions dominate. 

The numerical results were obtained by using the NLLBHA fortran code . 



The analytical and the numerical calculations for the cross-section in the non symmetrical 
Narrow- Wide configuration are in progress and will be presented elsewhere. 
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Appendix A 

Infinite momentum frame kinematics 



In this Appendix we will consider the kinematics we use to obtain the electron-positron 
and photon distributions. Due to the peculiar range of momenta and angles of the reaction, 
it is particularly convenient to use the Sudakov parametrization or infinite momentum frame 
kinematics. For the reaction 

e + (p 2 ) + e"(pi) -> e + (q 2 ) + e"(gi) + 7 (fc) (A.l) 

let us introduce the Sudakov decomposition: 

qi = a x p 2 + Pipi + qi, qi = a 2 p2 + P2P1 + q 2 
k = ap 2 + Ppi + k ± , (A.2) 

where p\^ are almost light-like four-vectors: 

q^Pi = qtto = 0, {qi? = -(^) 2 < 0, (A.3) 

2 2 
m _ m 

Pi = Pi P2, P2=P2 Pi, 

S S 

2 2 2 2 2 1 2 n ~2 ~2 m 

Pi = P2 = Qi = Q2 = m , k =0, p 1 = p 2 = — , 

s z 

s = 2p x p 2 = 2p 1 p 2 = 2p x p 2 = 2p 2 p x > m 2 , (A.4) 

where qf are Euclidean two-dimensional vectors in the center-of-mass reference frame. 

We consider the kinematical configuration when the photon is emitted in the direction 
close to the initial electron. We have the mass-shell conditions: 

q\ = sa 1 p 1 - (qi) 2 = m\ a, = ^ + m \ (A .5) 

spi 

{q' 2 f = sa 2 p 2 -{q^Y = m\ ft = + ™ 2 



k 2 = sap-(k ± ) 2 = 0, sa = 



sa 2 



P ' 

a 2 = 1, \f3 2 \ ~ I oil I ~ I ok j <C 1, Pi ~ (3 ~ I. 

The components along pi of the jets containing e~(qi) and ^(k) have a value of 0(1). The 
phase volume decomposition with d 4 gi = | daid(3id 2 qi is: 

d 3 q 1 d 3 q 2 d 3 fe .... , . . . „. 

d0 = 2 ; M L ^(pi+p 2 -qi-q 2 -k) (A.6) 

' d/3d/3i<J(l -0- pi)d 2 k ± d 2 qid 2 q 2 L 5 {2 \qi + q 2 L + k^). 



AsPP 
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The conservation law reads (we introduce a new four-momentum q of the exchanged pho- 
ton): 

Pi + q = qi + k, p 2 = Q2 + Q- (A. 7) 

The inverse propagators are (here and further we use f3\ = x): 

(pi - kf - m 2 = -^—d x , ( Pl + q) 2 -m 2 = — ^ -d, (A.8) 

1 — x x{± — X) 

g 2 = -(qi) 2 , d = m 2 (l - xf + (qi + q ± x) 2 , di = m 2 (l - x) 2 + (q{ + q 1 ) 2 . 
The matrix element reads 

M = -^v(p2)l^v(q 2 )u(qi)O u u(pi) 

°=1 7 tt- 2 2 e + e l — 2^1 ■ ( A - 9 ) 

(Pi — k,y — m z (pi + q) A — m z 

The following decomposition of the metric tensor g pu is used: 



9fiv = 9nv+ 1 + C( ) • (A.10) 

We use also the identity 

P2u(Qi)Ovu(pi) = u(q 1 )v p u(p 1 )e p (k). (A.ll) 
The generalized vertex v p has the form ||14|| : 

„^^ (1 -,)(I-i)-^(l-,)A(l-,). ,A, 2 ) 
The evaluation of the spin sum of the squared matrix element gives 

\v(Q2)Piv(p 2 )\ 2 = Tr P2P1P2P1 = 2s 2 , (A.13) 

spin 

The squared matrix element for the single photon radiation is given by 

R = -^Tr (p\+m)v p (p 1 + k-q + m)v p (A.14) 
= x[-2xm 2 (d - di) 2 + (q^) 2 (l + x 2 )dd x \—^ . 
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Finally we obtain that 



da 



e+e ^e+(e 7) 



= 2 a 3 d2q i dq i dx ( l 



X) 



\-2xm 2 (d - rfx) 2 + (qi ) 2 (1 + a; 2 )ddi]- (A.15) 



In the same way we may obtain the cross-section for the process of the double bremsstrahlung 
in the opposite directions: 



da ~~* e » e 7 
d 2 qrj L d 2 qr2"da;ida;2 



dV 



x 



x 



a\l + x 2 )(l + x 2 

- Xl )(l - x 2 ) J ((^) 2 ) 2 
(q ± ) 2 (l-x 1 ) 2 2xi m 2 (l-x 1 ) 2 (d 1 -d 2 ) 2 
did 2 

(^) 2 (1- 



(A.16) 



X2 



1 + xj d\d\ 
2x 2 m 2 (l-x 2 ) 2 (d 2 -d l ) 2 - 



d x d 2 l + x 2 d\d\ 



where xi, qi and x 2 , q 2 are the energy fractions and the components transverse to the beam 
axis of the scattered electron and positron, respectively; q^~ is the transverse two-dimensional 
momentum of the exchanged photon; 



[1 - x 1 ) 2 m 2 + (qi 



d x = (1 - xxfm 2 + (qi - q ± x 1 ) 2 , d 2 
di = (1 - x 2 ) 2 m 2 + (qi + q ± x 2 ) 2 , d 2 = (I - x 2 ) 2 m 2 + (qi + q L ) 2 



(A.17) 



Let us now discuss the restrictions on the d qi, d qi integration imposed by experimental 
conditions of the electron and positron tagging. We consider the emission of a hard photon 
along the electron direction. We will consider the symmetric case: 



0i<e e 

01<0e 



Q \ 

— <9 2 , 6 e =piq 1 , 
<6 2 , 6 S = p 2 q 2 . 



xe 
\Qi 



(A.18) 



Here 9\ and 9 2 are the minimal and maximal angles of aperture for the counters. It is 
convenient to introduce dimensionless quantities p = 2 /9i, z\ )2 = (g lj2 ) 2 /Q\ (Qi = 

The region in the 2:1,2:2 plane that gives the largest contribution to E is made by two 
narrow strips along the lines z\ = z 2 and z\ = x 2 z 2 . Therefore the leading logarithmic 
contribution will appear only in the cases where these lines cross the rectangle defined by 
x 2 < Zi < p 2 x 2 , 1 < z 2 < p 2 . Note that the line Z\ = x 2 z 2 , which corresponds to the 
emission of one hard photon along the momentum of the scattered electron, is the diagonal of 
the rectangle defined above. As for the line z 1 = z 2) which corresponds to the emission along 
the initial electron momentum, it crosses the rectangle only if x 2 p 2 > z 2 , xp > 1. This last 
condition defines the appearance of leading contributions to T, H . 
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For the contribution from the photon emission by the initial electron we have: 



p 2 x 2 p 2 



Fx = e(l-px)l^ f ^ziil-x 



z\ J (zx - xz 2 )(z 2 - z x ) 

1 X A 



p 2 x 2 p 2 



+ e(xp-i)/^ / 



dz 2 f dzx z 2 (l — x) 



z\ J (zx - xz 2 )(z 2 - Zx) 



1 2 

urpr Z2~r} 



(zx - xz 2 )(z 2 - Zx) J (zx - xz 2 )(zx - z 2 ) 

Z2+V 



d^i 2x<t r 2dz\ z 2 



2:2-17 « 2 -»7 

where we introduced the auxiliary parameter 77, a 2 77 1. Neglecting the terms of order 
rj we obtain: 



p 2 



Fr = f™{e(px-l)etf?-z)(L-^) (A.20) 
1 

+ e(x 2 P 2 -z)L 2 + e{z-x 2 P 2 )L 3 \, 

where Lj are given in eq. ( |3"TD and we used the identity 0(1 — px) + Q(px — 1)0(2 — x 2 p 2 ) = 

e(z-x 2 P 2 ). 

In the same way we obtain for the final electron emission: 



p 2 



dz ( , 2x 



zn £ -TT^ + £l )' (A ' 21) 

1 

The total contribution due to one hard photon emission in small-angle Bhabha scattering 
therefore reads: 



l-A 2 

X" = - / dx±±^-(F 1 + F 2 ). (A.22) 

7T J 1 — X 

Xc 
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Appendix B 



The contribution to E from the semi-collinear region of emission 
of two hard photons in the same direction 



An alternative way to use the quasi-real electron approximation is to compute the cross- 
section directly. To logarithmic accuracy we may restrict ourselves to considering only two 
regions i) the one with photon with momentum k\ emitted along the momentum direction of 
the initial electron inside a narrow cone with opening angle 6q <C 1, and ii) the region with 
the photon emitted along the scattered electron. Taking into account the identity of photons 
with the statistical factor 4 we obtain the cross-section: 



x 



r dV r dV " f 

J Tx((q^YYJ 7T J 



1-2A 



da; 



xix 2 x J 



(B.l) 



L-x- 

/ 



d^ 
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where 



f R ¥*t 2(0^0* [ d2fc ^ [l + a-^n^ + q-x!) 2 ] 

J 7T ^2)^1] n \ x ^i_ Xi y^p 1 k 1 )(2p 1 k 2 )(2q 1 k 2 

x[l + (l-x 2 f}[x 2 + {l-x 2 f] 1 
rri(l - x 2 ) 2 (2q 1 k 1 )(2p 1 k 2 )(2q 1 k 2 ) k^qj' 



fcillPi 



(B.2) 



(l-x 2 ) 2 (2q 

mai 
have 



Xi 

It is convenient to specify the kinematics: in the case of the emission 
>arallel to we have 



It is convenient to specify the kine: 
with momentum ki parallel to p 1 we 

2pih = 5i[(fe- 

Xi 

2qih = ^[xqi-(l- Xl )qt] 2 , 



Ql r/i-±\2 , 2 2 



a x 1 , 2pi/c 2 = — 

x 2 



of the collinear photon 
(B.3) 



in the case when the photon is emitted along q l we 

Q 2 



i,_L _1_ _i_ . 

K 2 ~ ~H2 "l 1 



: have 



2kiqi = ^L\a 2 x\ + (xki - qi)% ^2 = ^(H) 

X\X x 2 

1 Ql/ ± ±\2 i± ± J_ 1 — ^2 

'ih = (9i - zg 2 ) , fc 2 = -9 2 - 9i 

X 2 X 



(B.4) 



2gi 

= m 2 /Ql, and we introduced two-dimensional vectors k^, qf and q£ so 

^ 2 — QTirl r^7,± — sh 



where Q\ = e 2 6f, a 2 = m 2 /Q 2 , and we intrc 
that (q^) 2 = z u (q^) 2 = z 2 and q{q^ = <p. 
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The integration over d 2 k^ can be done with single logarithmic accuracy: 



Qi 



ir(2p 1 k 1 



feillPi 



= xiL, Q\ j 



7r(2q 1 k 1 ) 



feilWi 



—L. 

x 



(B.5) 



It is also necessary, here, to consider the kinematical restrictions on the integration variables 
4> and z\. When the photon is emitted within an angle 9q along the direction of the momentum 
of the initial electron, 9q represents the angular range to be filled by collinear kinematics events. 
We assign to the semi-collinear kinematics the events characterized by 



k 2 

X 2 



>0 , 



11) 



Qt 

X 



k 



2 

X-2 



>o , 



(B.6) 



where the region i) the photon with four-momentum k 2 escapes the narrow cone with opening 
angle 9 along the momentum direction of the initial electron. In the region ii) the same 
happens for the final electron. 

We can rewrite the conditions above in terms of the variables Z\ and <fi as follows: 



i) 
ii) 

Hi) 



1 > COS0 > -1 + 

1 > COS0 > —1, 

1 > COS0 > -1 + 
XJ1~2~ 



A 2 -(v^r-y^) 2 

2^/Z\Z 2 

\zi - z 2 \ > 2 v /i 2 "A, 



y/zi- \fz2\ < A, 



2y/ziz 2 — 



■XI 



iv) 



1 — X\ 

1 > COS0 > —1, 



< A— 
1 

w - 



X 



Xi 

X 



X . X 

— —Z 2 \ > 2X v /z 2 -7- r^, 



where A = x 2 9 /9i. In our calculation we take the parameter A < 1. Indeed, the restrictions 
on #o for collinear kinematics calculations are e9 > m or O > 1CT 5 at LEP energies. On the 
other hand the experimental conditions on 9 1 are 9\ > 1CT 2 . Therefore we can take A C 1 
within our accuracy. 

Analogous considerations can be made for the case when a photon with momentum k\ 
is emitted along the direction of the final electron. In regions ii) and iv) we may do the 
integration over the azimuthal angle: 



2vr 



dc/> 







2n(2p 1 k 2 )(2q 1 k 2 ) 
d0 



x 2 xQ x 4 



lo 2ii{2pik 2 ){2q l k 2 ) 



feillPi 
feilWi 



x(l — Xi 



1 - x 1 )z 1 - xz 2 l\z 2 - Zi\ 



x 2 x 3 (l - x 2 ) 2 Qi 
Zi - z 2 x 2 /(l - x 2 ) 



\x 2 z 2 — (1 — xi) 2 zA 



1 - x 2 



\zi-Z2jT^\ \zi-x 2 z 2 \ 



,(B.7) 
.(B.8) 
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The integration of regions i), Hi) has the form 



1 



7T J \Z\— Z 2 \ 



d0 

2tc(z 1 + z 2 + 2^/ziZ 2 cos0) 
\zi - z 2 \ 



dz 



arctan 



tan - 



(B.9) 



where 



) = arccos 



2^/z\Z 2 



The result reads 



J = 2 In 2. 



(B.10) 



(B.11) 



We give here the complete contribution of the semi-collinear region: 



da. 



HH 

coll 



q 2 £ 
4tt 2 



1-2A 1-x-A 



da; 



dxidx2$(l — x — X\ — x 2 ) 



x 



£1X2(1 - ^i) 2 

2 

/ ^{ln ^[1 + 6(pV - *) + 26(p 2 (l - Xl f - *)] 
1 

(z — x 2 )(p 2 x 2 — z) 



[l + (l-£ 1 ) 2 ][x 2 + (l-a; 1 ) 2 ] 



+ 9(pV-z)m 



x 2 [z — x(l — Xi))(p 2 x(l — Xi) — z) 



+ Q{z-p 2 (l-x l f 



In 



(z — p 2 (l — xi)x)(z — (1 — Xi) 2 ) 
(p 2 (l — xi) 2 — z)(z — x(l — Xi)) 



+ In 



(p 2 (l-x 1 )-z)(z-(l-x 1 ) 2 ) 
{p 2 {1 _ Xl y_ z){z _ {1 _ Xl)) 



+ 0(z- p 2 x 2 )\n 



z — p 2 x(l — Xi)(z — x 2 ) 
(p 2 x 2 — z)(z — x(l — Xx)) 



+ Q{p 2 {l- Xl ) 2 -z) 



In 



(z-jl-x^jp^l-x^-z) 

( p 2 a .( 1 _ Xl) - Z )( Z - X (l - Xl ))(l _ Xl )2 



+ In 



( Z -{1- Xl f)(p\l- Xl f- Z ) 



+ In 



^-1)(P 2 -^) 



z _ ( l_ Xl))(p 2 (l _ ;ri) _ z) 



(p 2 (l-x 1 )-z)(z-(l-x 1 ))(l-x 1 ) 2 . 

To see the cancellation of the auxiliary parameter #o/#i we give here the relevant part of the 
contribution for the collinear region : 



1-2A 1-x-A 



a 

4^2 



da; 



dxi da; 2 5(l — x — xi — x 2 ) 



[l + (l-£l) 2 ][£ 2 + (l-£l) 



X 



l+ 1 -e(p 2 x 2 -z) + Q(p 2 (l-x 1 ) 2 -z) 



+ ... . 
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We see from the above expression that the dependence on 9q/9i disappears in the sum of the 
contributions for the collinear and semi-collinear regions. The total sum is given by Eq. fl56|). 
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Appendix C 

Virtual corrections to single photon emission cross— section 

The cross-section for single hard photon bremsstrahlung containing virtual and real soft 
photon corrections may be written as follows: 

d^+v) = ^VdV R= lim (C1) 

2n 2 x(l -x){q^Y ' (2 PlP2 Hoo (2p 1 p 2 ) 2 y ! 

We define the Compton tensor with a heavy photon as [17]: 

spins 

where M p is the matrix elements of the Compton scattering process 

7*(<Z) + + e(p2), (C.3) 

i 2 n 2 2 2 

fc = 0, Pl =p 2 = m , 

p is the polarization index of the heavy photon (q 2 ^ 0). The Bose symmetry and the gauge- 
invariance requirements provide the following general form of K pa : 

K pa = ^(P PCT + P; p ), (C.4) 
Ppa = 9pa(B g + ^-Tg) + plppla(B U + ^ T ll) + P2pP2a(B u + ^-T U ) 
+ PlpP2a[Bl2 + ^-T 12 ) + P2pPla{B 12 + ^~T 12 ) , 



where 



Z7T Z7T 

g P P P <7 = 0, q a P pa = 



HpHcr ~ VII ~ V21 tri r\ 

9 P a = 9 P a - , Pip = Pip - Qp—^- , P2 P = P2 P - q P ~^ ■ (G.5) 

Coefficients Bi and 7] before different tensor structures depend on the kinematical invariants 

s = (p 2 + k) 2 - m 2 , t 2 = (pi - k) 2 - m 2 , u = (p 1 -p 2 ) 2 , s + t + u = q 2 . (C.6) 
In the Born approximation we have 

B g = j t [{s + u) 2 + {t + u) 2 ]-2m 2 q 2 [^ + ^y B 12 = 0, 

— B i = Bi(s,t,u), Bi = Bi(t,s,u), i = g, 12,11. (C.7) 



4g 2 8m 2 
-on = — - 
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Contributions in the one-loop approximation are: 



12 



11 



+ 



+ 



r = 



rB g + 



1 (u 
s\s 



f + \ i 2ub + s 2 ))G + u( q ---\l q -l u ) 



st a 



(C.8) 



c /3m 



sV b 



-l)(l q -h) + 



u 



s m 
2st + ~W 



m 2 s 



'-.(4b -u) -2s )h 



2u + AbG - 2q 2 (l 2 - V u ) + 2( b + — )F - bn 



-\-Au - s)G + ^ (uq 2 - st)G - 2q 2 ' ""' 
st s z t z 



f (t <-> s + iO) j, 

(/g Z M ) 



■uo 2m — s + t 
st a 



- « 

a 



cs) (^(Z, - Z u ) - l) + ^ (2c + t) (s - (l g - l s ) 



Cp 1 QiUjS 

±— (2m - s)(L - Z t ) + 8m + 3t - s + 

bs c 



+ m 



t 2c 

^(ut + 2m 2 b)l t + 2[ - + — )l s - 



u + 2b u + Ac 2c ~ 



t 



b „ c ~ s+6 m 

+ -.R+-.RH n n 

u u t s . 



u 



st 



rB n + -{-q 2 [l + - )G - q 2 [2 + - )G + 2q 2 [- + — )(l q -l u 



+ —N + Q 

Sa 



2u 



st a 

2 



A 2 2 z,2 2 

+ -(u 2 - bt) (?-(l q - Q - l) + ±— (2c + - Z s ) + ^-(2u- s)(l q - l t ) 



4m - 2g 2 + i Ym 2 

c 



, 41t - 5)+2 ((^)!-f), 



2m 



2 + 



m 



8t 
s 
b 2 



—N -Q- —R -R + -n (2c 2 - u 2 )h 



S a 



A 



t cs 

7T 2 9 



4(Z u -l)ln--Z 2 + 3Z (Z + , 2 . 



where the following notation is used: 

i = t + m 2 , s = s + m 2 , a — s + t, b = s + u, 

9 2 

m m 

c = t + u, s a = s a, t a — t a, 

u u 

G = (Z,-Z u )(Z, + i„-2Z t )-2Li 2 (l)-2Li 2 ^l-^+2Li 2 ^l-^, 
F = _Li 2 (l + -V)+Li 2 (l), N = 2l t (l q — l u ) + 2Li 2 (l — — ^ — 2F, 



(C.9) 



u 
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^ Ac ( m 2 \ ~ A m 2 ^ m 2 ft, \ 

q = -^--y-u-F, „ = T ( |ll -i), 

it = lnf — /g = ln(-^-], l u = \n( — ^-V f s = In f . 

In order to avoid the dependence of the Compton tensor on the photon mass A we need 
to take into account the emission of an additiona soft photon with energy uj < As. The 
corresponding contribution, which has to be added to the quantity r, in the case of small 
angle Bhabha scattering reads [15]: 



771 IT 

p = 2{l u - 1)(2 In — + 2 In A - lux) + l 2 u - In 2 x 

A 3 



(CIO) 



Note that in the case of small-angle Bhabha scattering the contributions of the g pa tensor 
structure is suppressed by factor 6 2 relatively to terms containing pi P pj a . 

We need only limited values of T ik in the cases of s 1 <C \ti\ and \ti\ <C s 1 at fixed q 2 and 
u\ = —2piqi. In the case of small Si we have s± = s — [m 2 (l — x) 2 + (q^x + qf) 2 ]/[x(l — x)] 
(we omit in the remaining part of this Appendix the subscript 1 in the notation of invariants 
of the Compton subprocess). 

Taking into account that, at small s, q 2 = — (q^) 2 , t = —(1 — :r)(g^) 2 and u = — (q^Yx, 
we derive the following expressions for h = p + r and T in this limit: 



h s<t \ = 2(L-l + lnx)(21nA-ln:r) + 3L-ln 2 



x — — 



9 



2' 



s<|*| 



s(l - x) 
l + 2x + x 2 



4(1 +x 2 

16m 2 



(a 1 ) 2 

lux In - Li 2 (l - x) 



InxL. 



In the case of small \t\ we have: 



(C.11) 



h ltKs = 2(L-l-lnx)(21nA-lnx) + 3L-ln 2 x 



2' 
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|t|<s 



2X f . n 



(q r ) 2 1 
lnx In In x — Liofl — x) 

-t 2 y . 



1 - 2x + x 2 } + 



16m 2 x 2 
~1F 



lax L. 



(C.12) 



One can be convinced that all terms reinforced by large logarithms (L or L 2 ) arise only from 
these regions. We omit here only terms of order (a/n) 2 . 
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The further integration is straightforward. We show here the most important moments. 
The contribution of the h containing terms gives (in close analogy with the Born contribution): 



^H(S+V) 



2 ir' 



1-A 



/W d 4^{( i+e (" v -)> 

1 x c 



(C.13) 



3 1 15 

x L(21n A - lux + -) + (2 In A - lux) (In a: - 2) - -In 2 x - — 

3 1 
+ (21nA-lnx + -) k(x,z) - 21nx(21n A - lna;)6(p 2 :r 2 - z)\. 

To obtain the contributions from T we consider at first new types of integrals: 

i s {t} = Ql 
i s {t} = Ql J 



fd 2 qi ln (qiy 



n(qir J 7rs{t} s{-t} ' 



(C.14) 



d q.2 f d <ii m2 



n(qiyj irs 2 {t 2 } 



n(qir J 7rs{t} ' 

Denoting <r 2 (l — x) 2 + (q 2 x — qj;) 2 /Q 2 as a + 6cos0 and using the expressions 

2tt 



1 



d0 



1 



27r J a + bcos< 



2tt 



1 /" ^hi(a + 6cos(/>) 
+ 6 cos <fi 



27r J a 



Va 2 - 6 2 ' 



In 



(C.15) 



2 (a 2 - b 2 



Va 2 - b 2 a + Va 2 - b 2 



with a 2 — b 2 = (zi — x 2 z 2 ) 2 + 4<r 2 (l — x) 2 x 2 z 2 , o 2 = m 2 /Q\ and z^ 2 = (qi, 2 ) 2 /Ql we derive 
that 



d 

I s = (1 - x)x J j dzi 



(C.16) 



x 



lnQz 2 (l - x)x 3 ) - ln[(zi - z 2 x 2 ) 2 + 4aV(l - xfz 2 \ 
(z 1 - x 2 z 2 ) 2 + 4cr 2 x 2 (l - rr) 2 ^ 



Since we are evaluating with logarithmic accuracy, we may consider the contribution of the 
region \z\ — x 2 z 2 \ < rj, a 2 <C rj <C 1 when integrating over zi. The result reads 



= x(l — x) J 



dz 2 

T 2 " 



1 T X 

-L + In 

2 1 -x 



(C.17) 
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The rest of the integrals can be calculated in the same way, and we have: 



p 2 x 2 



-a-x) 



-L 



l T , 1 

-L + In 

2 l-x 



dz? 



it = -(l-x) / -£L, (C.U 



i 



x{ 



l 2 l 



p 2 x 2 p 2 X 2 
d£2 2 

zi 



dz 2 

z 2 



Using ( C.17 ) and (|C18 ) we may represent the final result for the contribution to '£ H ( S+V ) due 
to the T term as 



1 , a , o P r dz _ 1 /" A , 1 + x 2 ( ( 1 



2V J z 



L dx 



l-x LV2 



-L + ln 



1 — x 



\nx 



Xc 

„2 



x +2x — 1 2x lnx 



+ C2 - Li 2 (x) + — ^ - - 6(xV - z) 



1 



4(1 + a; 2 ) 1 + x 2 

1 + 2x — x 2 2x Inx 



2 l-x 



(C.19) 
lnx 



— - In 2 x — Liofl — x) — 

2 v ; 4(1 + x 2 ) 1 + x 2 

The total contribution to J2 H ( S+V ^ (one-side hard photon emission with virtual and soft 
photon corrections) is the sum of ( C.13|) and ( P.19 ): 



j^H(S+V) __ j^H{S+V) _|_ j^H(S+V) 



(C.20) 



This quantity is given in Eq. 
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Appendix D 

Leading logarithmic contribution to £ 77 

Here we show that the main logarithmic terms can be summed according to the renormal- 
ization group. The sum, as given in the text, may be written as: 



(D.l) 



x (i + e(*y-z)) + - 



1 /-i-A 1 + X 



2 In A - lnx + -)(1 + 9(>y - z)) 
1 — x 2 



1-2A 



+ 



dx 



x 2 , 1- x- A 1 



1 — x 



In 



A 



1 + x)lnx — 1 + x 



[l+3Q(xY-z)] 



+ 



dx 



1+x 2 

1 — X 
-,2 



In 



A(y/Z- PX) 
Z X P , ■> , 2 2,1 

B(2; — x p ) >. 



sfx + x - 1 



- i(l + x)ln— + 

2 z p 

One can see that the dependence on the parameter A disappears in the expression above. 

We will now show that eq. ( |L). 1| ) is equivalent to eq. (|59|). Let us transform eq. using 
the substitution 



Q{t 2 p 2 - z ) = 1(1 + 6(a; 2 p 2 - z)) + ^9(z - x 2 p 2 ) - &{z - t 2 p 2 ) 



(D.2) 



and changing the order of integration in the last term: 



J p z 
dt f dz 



dz 



dt 



dzQ{z - p 2 x 2 ) 



dt. 



(D.3) 



pH 2 



By evaluating the integral over t, and using the explicit expressions for the splitting func- 
tions one can verify the coincidence of Eqs. ( p.l|) and (|59|) . In an analogous way one can prove 
the validity of the representation ( |61~D for S 7 . 

Using the representation in Eq. ( pOD for the function one can see that the above 
expression is equivalent to Eq. (p.l|) . In an analogous way one can prove the validity of 
representation in Eq. floT| ) for E 7 . 
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Appendix E 

Cancellation of the A dependence in the nonleading contributions to £ 27 

Let us consider the singular nonleading terms in E 27 in the limiting case A — > 0. Dropping 
the common factor (a/n) 2 C J dz/z 2 , we give below the various contributions separately. 

Let us consider first X 77 . The contributions from the soft photon radiation and virtual 
corrections are: 



.vv+vs+ss^ = i n A(-7-41nA)(l,p 2 ), (E.l) 

where we denote by (a, b) the limits of the integration over z: (a, b) = Q(z — a)<d(b — z) . 
The contribution due to the virtual corrections to the single hard photon emission gives^ 



{E h(s+v) )a = li ri A|(l,p 2 )(161nA + 14)+4|dx^^[(l,p 2 )-(l,p 2 x 2 )] (E.2) 

Xc 

1 1 
+ 4[-2 ln(l - x c ) - 2 ln(l - x c ) + J dx (1 + x) + J dx (1 + x)} 

X c X c 

1 1 2 1 2 

+ 2 J dx T ^[(l,p 2 )-(l,pV)]lna; + 2 J dx k(x, z)}, 



where x c = max(x c , 1/p) and the quantity k(x, z) is defined in Eq. (|30|) . The singular part in 
the contribution in Eq. (EB) due to double hard photon bremsstrahlung reads: 



2 

(Z HH ) A = lnA{|dx T ±^[-(l,p 2 )L 1 -(l,pV)L 2 -((l,p 2 )-(l,pV))L 3 ] (E.3) 
1 1 

- J dx (3 + x) - J dx (3 + x) - 4 In A + 4 ln(l - x c ) + 4 ln(l - x c ) 

- /d^i±^[(l, P 2 )-(l,pV)]}. 

Xc 

It is possible to verify the cancellation: 

^vv+vs+ss )a + ( E h(s+v ))a + (s ^h )a = (E4) 

9 In our contribution to the Tenessee-94 workshop (see [2]) there are some misprints, which are corrected 
here 
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The corresponding contributions to are: 



(Ef#) A = lnA(-14-81nA)(l,p 2 ), 

1 2 

(Sf +V + Sf fy ) A = lnA^/dx^- A;(a;, 2) + (1, p 2 ) [16 In A + 14 - 8 ln(l 



1 



- 81n(l -x c ) +4jdx (1 + x) +4y dx (1 + x)] 

+ 4/dxi^[(l,p 2 )-(l, P V)]}, 

(Eg) a = lnA{-8(l,p 2 )[lnA-ln(l-:rJ-ln(l-iJ] 

_ g/ <to (l.^)-(l.^)_ 2 / djt( Xtt) h±* 
J l-x J v 1 1 

1 

2\ i /i „2 2\ 



- 4 j /d:r(l + :r)[(l ) p 2 ) + (l ) pV)]j. 

Rearranging the last term in (£#)a as 

i i 
-4 y dx (1 +x)[(l,p 2 ) + (l,pV)] = -4jdx (l + x)(l,p 2 ) 

X c X c 

l l 
-4 / dx (1 + x)(l,p 2 ) +4 / dx (1 +*)[(!, p 2 ) - (l,pV)], 



Xq Xq 



we can see again the cancellation of the A-dependence in the sum: 

(EfSU + (Sf +V + E| +v )a + (SI)a = 0. 
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Appendix F 

Relevant integrals for collinear pair production 



We give here a list of the relevant integrals, calculated within the logarithmic accuracy, for 
the collinear kinematical region of hard pair production. 

We use the definitions in Eq. fl6T|) and we imply, in the left-hand side of the expressions 
below, the general operation: 



2tt_ 



«...»./-*/-/£(...) 



(F.l) 



with the conditions z = (e9 /m) 2 3> 1, L — lnz 3> 1. The details of the calculations can 
be found in the Appendix of Ref . |B| . The results are: 



,(x 2 D + (1 -x 2 )A" 



V DC J ' 

(l-x)\l-x 2 f 



x 2 ) 



L + 2 In 



XiX 2 



X 



(x 2 A x - x\A 2 
V AD 



XX\X2 
2 



(— >=- 

DC Xix 2 (l — x 2 ) 2 



[-Lq + In ■ 



} 



[l-xf 



L Q + 2 In 



x x x 2 

X 



'1-x) 



XX\X 2 



X 

4(1 -x) 

X 



x\A 2 - x 2 Ai _ {x\- x 2 )L 



AD 2 ' xx\x 2 (l — x) ' 

—L n r l _ , . X\X 2 . 



AD' 


0:1X2(1 — x) 2 


1 ) 

AC 


= -%L + 21n 

X\X^ 


A 


v ^2-^0 


C 2 D'< 


!/ zi(i-z 2 ) 4 ' 


A 


)= " Lo [- 


CD 2 


7 X!(1-X 2 ) 2L 2 


C 

A 7-19 


»=.. - L ° [i 



(-L) = ~ L ° 

1 



X 



X\X 2 



+ 21n 



( C2D ) 

xx 2 



-L 



Xi(l - x 2 ) 3 ' 



X 



<— > 

h m 

a; 



(1 - x)(l - x 2 ) 
-x 2 L 



-0 



XXiX 2 



X) 



X 



\ + L 
-Lq\ 



X 2 X — X\ 



XX\X 2 {\ — x 2 ) 2 
' X\ — x 2 



X\X 2 {1 — x) 2 xx 2 (l — x) 



(F.2) 
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